Isospin Dependence in the Odd-Even Staggering of Nuclear Binding Energies* 
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The FRS-ESR facility at GSI provides unique conditions for precision measurements of large ar- 
eas on the nuclear mass surface in a single experiment. Values for masses of 604 neutron-deficient 
nuclides (30<Z<92) were obtained with a typical uncertainty of 30 fiu. The masses of 114 nuclides 
were determined for the first time. The odd-even staggering (OES) of nuclear masses was system- 
atically investigated for isotopic chains between the proton shell closures at Z=50 and Z=82. The 
results were compared with predictions of modern nuclear models. The comparison revealed that the 
measured trend of OES is not reproduced by the theories fitted to masses only. The spectral pairing 
gaps extracted from models adjusted to both masses and density related observables of nuclei agree 
better with the experimental data. 



Significant progress has been achieved over the last 
years in constructing self-consistent mass models 0, . 
These models aim to reliably describe the properties of 
nuclei far off the valley of /3-stability, where the exper- 
imental information is scarce or even not available yet. 
For instance, in modelling the astrophysical r-process of 
nuclear synthesis one needs precise knowledge of masses 
and half-lives of very exotic nuclei and one has to rely 
on theoretical predictions since most of the nuclides in- 
volved have not even been produced in the laboratory yet. 
The predictions for these nuclides dramatically deviate 
for the different models [2j ■ Thus new experimental data 
on exotic nuclei and consequently better understanding 
of nuclear structure away from the valley of /3-stability is 
essential for further theoretical development. 

Odd-even staggering of nuclear binding energies (OES) 
was detected in the early days of nuclear physics [3j and 
was explained by the presence ofpairing correlations be- 
tween nucleons in the nucleus |4j. Pairing contributes 
only little to the total nuclear binding energy but its in- 
fluence on the nuclear structure is significant. 

The common way to extract experimental information 
about the pairing correlations is to measure the value 
of the OES which approximates the pairing-gap energy 
(A) in the standard Bardeen-Cooper-Schrieffer (BCS) 
theory ||. The latter quantity is connected with the 



strength of the pairing interaction (G): 
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(e„ - A) + A 2 
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where e„ is the single-particle energy and A is the chemi- 
cal potential. The summation goes over all single-particle 
levels v below and above the Fermi energy. In order to 
evaluate this sum in local pairing functionals a (smooth) 
cut-off in energy is usually implemented. 

Neutron (A„) and proton (A p ) pairing gaps are usually 
determined from finite-difference equations of measured 
masses [(j, e.g. by the five-point formulae: 

A,^ =-- [M(Z, N + 2) - iM(Z, N + 1) + 
8 

6M(Z, N) - 4M(Z, N - 1) + M(Z, N-2)], (2) 

= -i [M(Z + 2, N) - 4M(Z + 1, N) + 
GM(Z, N) - AM{Z - 1,N) + M(Z - 2, N)], (3) 

where M(Z,N) is the mass of an atom with Z protons and 
N neutrons. 

The well-known parametrization A~ Ylj\fA MeV 
(A = N + Z) provides the average trend for nuclei close 
to stability. A dependence of the pairing strength on 
the neutron excess was suggested in @. It was later 
observed from the mass determination of exotic Dy-Hg 
isotopes that A p and possibly A„ increase towards the 
proton drip-line 0. 
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In this letter we present new results on the OES ob- 
tained from our high-precision mass measurements com- 
pared with predictions of modern nuclear theories. 

The experiment for direct mass measurements was 
performed at the FRS-ESR facility as continuation of 
a successful scientific program addressing basic nuclear 
properties of stored exotic nuclides 0, E3 EH- Exotic 
nuclei were produced by projectile fragmentation of a 
(600-900) MeV/u 209 Bi primary beam in (4-8) g/cm 2 
beryllium targets placed at the entrance of the frag- 
ment separator (FRS) [l2j. The fragments were spa- 
tially separated in-flight and injected into the cooler- 
storage ring ESR E3l- In the ESR, the velocity spread 
of the stored fragments was reduced by electron cooling 
to Sv/v w 5 • 10~ 7 . This condition provides an unam- 
biguous relation between the revolution frequencies of the 
ions and their mass-to-charge ratios which is the basis for 
Schottky Mass Spectrometry (SMS). The time required 
for the electron cooling was about 10 s, which constrains 
the range of nuclides that might be investigated by this 
method. The SMS has reached the ultimate sensitivity by 
recording single ions stored in the ESR leading to a mass 
resolving power of more than 2-10 6 (FWHM) [3 Ell- 
in this new experiment with SMS, 582 different nu- 
clides were observed in the frequency spectra. From this 
set of nuclei, 117 were used for calibration. The achieved 
mass accuracy was typically 30 \±u which represents an 
improvement by a factor of three compared to our for- 
mer experiments E3- I R addition, the masses of 139 
nuclides were determined indirectly by means of known 
decay energies (a, /3, or proton emission). The masses of 
114 nuclides were obtained for the first time EH- The 
measured masses cover a large area of neutron-deficient 
nuclides from krypton to uranium. All directly measured 
values (see EH) have been included in the latest Atomic 
Mass Evaluation E3- 

The achieved experimental mass accuracy has allowed 
us to perform new systematical studies on nuclear pair- 
ing. The new data were combined with the data of 
Ref. E3 an d precise values of OES for all even-Z iso- 
topic chains in the region between the Z=50 and Z=82 
closed shells were extracted. Only even-even nuclei were 
considered. The results obtained show that indeed the 
values of OES for both the protons and the neutrons 
increase towards the proton drip-line, thus confirming 
earlier observations for a small number of nuclides [8j. 
Moreover, this is a general trend for all even-Z isotopic 
chains from tin to lead. The tin, tellurium, mercury, 
and lead isotopes were not considered for protons since 
the closed shells at Z=50 and Z=82 have strong influ- 
ence. Similarly, the nuclides with N=80, 82, 84, 124, 
126, and 128 were excluded for neutrons. It is necessary 
to note that no such general trend of the OES was ob- 
served for isotopes below tin. For isotopes above lead 
the experimental information is still too scarce to draw a 
definite conclusion. The large number of newly obtained 




FIG. 1: Comparison of the proton and neutron pairing-gap 
energies for even-even hafnium (upper panel) and tungsten 
(lower panel) isotopes derived with 2nd-order mass differ- 
ences of experimental masses and from the predictions of the 
original FRLDM mass model [T^ and with newly readjusted 
pairingstrengths. The experimental values are taken from 
Refs. El El. 



OES values allow us to perform quantitative comparisons 
with calculations. We will first compare the results with 
a macroscopic-microscopic model and then continue with 
up-to-date microscopic models. 

The pairing-gap energies were extracted from the 
masses calculated with the Finite-Range Liquid-Drop 
Model (FRLDM) El- In this model, the single-particle 
potential generated by the Yukawa interaction was used 
for the microscopic part. To obtain a maximal num- 
ber of extracted OES values, the calculations were done 
using 2nd-order mass differences (three-point formulae 
d). The comparison with the experiment is shown for 
the isotopic chains of hafnium and tungsten in Fig. ^ 
It is clearly seen that the experimental isospin depen- 
dence of pairing- gap energies is not reproduced by the 
original FRLDM [19|. To improve this description, the 
BCS pairing part of the model has been adjusted to 
the new experimental data. Different from Ref. E3j 
a single-particle spectrum was generated with the de- 
formed Woods-Saxon potential. The pairing strength 
G was parameterized with 2 constants for protons (p) 
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FIG. 2: Comparison of the proton (A p ) and neutron (An) 
pairing-gap energies for even-even hafnium isotopes derived 
from fourth-order differences of experimental and of calcu- 
lated masses. Hartree-Fock plus pairing treated with the BCS 
formalism (MSk7) and three HF-Bogoliubov models with dif- 
ferent cut-off parameterizations (BSkl and BSk2) and density 
dependent pairing (BSk3) were used. The experimental val- 
ues are taken from this work (full symbols) and from Ref. [Ts| 
otherwise. 



FIG. 3: Comparison of the proton ({A) p (uv)) and neutron 
({A) n (uv)) pairing-gap energies for even-even hafnium nu- 
clides calculated with several RMF and Skyrme-Hartree-Fock 
models (see text). These models were adjusted not only to 
the nuclear binding energies but also to form-factor related 
observables. The experimental points were derived from mea- 
sured masses using Eqs. I2I3H . The full symbols represent our 
new experimental values, the others are from Ref. [l^l . 



and neutrons (n) G p{n) = g 0p (n)/A + gi p ( n )(N - Z)/A 2 
0, LI- The number of levels taken into account (see 
Eq. 2]) was equal to N for neutrons and Z for pro- 
tons. All experimental OES values for even-even nu- 
clides between Z=50 and Z=82 were used for the ad- 
justment. The best agreement between the values of 
OES derived from the calculated masses and experimen- 
tal values was achieved withoop = 9on = 20.80 MeV and 
9ip = -gin = 22.40 MeV [3. Note that the constants 
(ffop, gon) and (gi P , gin) converge to the same values, 
which is quite remarkable since this was not a constraint 
demanded in the analysis. The results obtained are la- 
belled in Fig. ^ as readjusted FRLDM. The difference be- 
tween the original FRLDM of Ref. and the readjusted 
model is obvious. With this new pairing description the 
o rr ns value for the prediction of nuclear binding energies 
between Z=50 and Z=82 closed shells has improved by 
about 25 %. 



Going to microscopic models, pairing-gap energies 
were calculated from predictions of nuclear mass cal- 
culations of several self-consistent mass models, which 
are Skyrme-Hartree-Fock (SHF) calculations plus pair- 
ing treated in the BCS (HF+BCS) formalism (MSk7 
force) |2PJ and three models where pairing is treated with 
the Bogoliubov (HFB) approach with different pairing 
cut-off parameterizations (BSkl, BSk2) 12 lL |2'J and in- 
cluding density-dependent pairing (BSk3) [23j. The re- 
sults calculated with 4th-order mass differences (Eas. l2l3|l 
are presented for even hafnium isotopes in Fig. [21 It is 
clearly seen that the general trend is not reproduced by 
any of the models. A comparison with four other recent 
HFB models where the implementation of different effec- 
tive masses, with and without density-dependent pair- 
ing (BSk4-BSk7) showed similar results as those of 
Fig.H 

Another branch of self-consisistent mean-field mod- 
els are relativistic mean-field (RMF) models employing 
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finite-range meson fields (FR) or point couplings (PC). 
For this study, we employed three of the best param- 
eterizations available, namely NL-Z2 Eit and NL3 j2|| 
for the finite range variant, and PC-F1 [27| for the point- 
coupling model. In contrast to the mass models described 
above, these RMF forces are adjusted to both energy and 
form-factor related observables (e.g. rms radii, diffrac- 
tion radii, surface thicknesses, etc.), and are meant to 
describe both kinds of observables. Furthermore, we per- 
formed calculations with the SHF forces SkI3 and SkI4 
|28j with extended spin-orbit terms, which - similar to 
the RMF forces - have been adjusted to both masses 
and density related observables. In both the RMF and 
the SHF models, we employ BCS pairing with a density- 
independent (5-force. As a first approach, the pairing 
gaps have been estimated from the single-particle spec- 
trum with uu-weighted single-particle gaps |2^| (v 2 are 
the occupation probabilities), which circumvents the un- 
certainties related to the calculation of odd-even systems. 
These quantities constitute a measure of the pairing con- 
tribution to the OES. As discussed in Ref. j^, these 
results need to be carefully interpreted due to polariza- 
tion effects and the non-pairing-type contributions to the 
OES. 

It is striking that the RMF and SHF calculations in 
Fig. |31 give very similar results: the general trend of the 
rising pairing-gap energies is reproduced. However, some 
local discrepancies are observed, as e.g. close to T z =5, 
which can be related to the N=82 closed shell. Although 
the models in Fig. |3 in general have much higher pre- 
dictive power for the nuclear masses Q the difference 
in the description of the experimental OES data is obvi- 
ous. This result has not been expected and demonstrates 
the need for a better understanding of both the roles of 
the various observables and adjustment protocols as well 
as the procedure of calculating OES within these frame- 
works. 

All models tested in this letter take into account the 
nuclear deformation which is essential here since most 
of the nuclei investigated are deformed. Moreover, the 
observed general trend of OES has the same slope and 
magnitude for nearly all isotopic chains from xenon to 
platinum. Use of the 2nd, 3rd, or 4th-order mass dif- 
ferences to disentangle the mean field contributions to 
the OES and pairing-gap energies has been intensively 
discussed in the literature Our over- 

all conclusions, as checked in different analyses, are not 
changed if 3, 4, or 5-mass formulae are used. Since the 
mass number changes within an isotopic chain, volume 
effects might contribute to the observed isospin depen- 
dence of the OES. 

With the new data available, it became possible to 
examine the OES predictions of different theories. The 
new results are helpful for a better description of the 



pairing in exotic nuclei, which is mandatory for a reliable 
theory. 

An important future aspect is whether the observed 
trend of the pairing-gap energies persists for nuclides with 
even greater neutron excess. A recent experiment has 
been performed at the FRS-ESR to measure masses in 
the Yb-Pb region on the neutron-rich side of the chart 
of nuclides. Measurements of very exotic neutron-rich 
nuclides, which cannot be produced with the present fa- 
cility, are foreseen within the FAIR project |35| . 
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